Field-Induced Freezing of a Quantum Spin Liquid on the Kagome Lattice 
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We report ^^O NMR measurements in the S = 1/2 (Cu^'*') kagome antiferromagnet Herbert- 
smithite ZnCu3(OH)6Cl2 down to 45 mK in magnetic fields ranging from 2 T to 12 T. Wliile 
Herbertsmithite displays a gapless spin-liquid behavior in zero field, we uncover an instability to- 
ward a spin-solid phase at sub-kelvin temperature induced by an applied magnetic field. The latter 
phase shows largely suppressed moments < 0.1/iB and gapped excitations. The H — T phase di- 
agram suggests the existence of a quantum critical point at the small but finite magnetic field 
/io-ffc = 1.55(25) T. We discuss this finding in light of the perturbative Dzyaloshinskii-Moriya in- 
teraction which was theoretically proposed to sustain a quantum critical regime for the quantum 
kagome Heisenberg antiferromagnet model. 



Quantum spin liquids (QSL) are appealing states of 
matter which do not break any symmetry of the spin 
Hamiltonian. While QSL behaviors are well estab- 
lished for one-dimensional spin systems, their existence 
in higher dimensions remains questionable. The quantum 
kagome Heisenberg antiferromagnet (QKHA) is consid- 
ered as the best candidate to stabilize such a QSL phase 
in two dimensions [ll-0|- The specificity of this model 
relies on the unique combination of strong quantum fluc- 
tuations enhanced by low spins S ~ 1/2 and high geo- 
metrical frustration of the lattice of corner-sharing tri- 
angles. Early numerical studies have shown, despite fi- 
nite size limitations, that the ground state of the QKHA 
model does not support any simple order parameter 0- 
0] and have evidenced an original excitation spectrum 
with dense sets of low energy excitations in all spin sec- 
tors d, Q. A recent state-of-the-art calculation [7] fa- 
vors a QSL ground state with a small gap much like 
the resonating- valence-bond state ^, 9]. Other recent 
proposals with similar ground-state energies encompass 
large unit-cell valence-bond-crystals [l^l as well as gap- 
less 'critical' spin liquids with algebraic spin correlations 
111 Il2l | . Clearly, the exact ground state of QKHA re- 
mains a challenging and highly debated issue 0,113. In 
this respect, studying the effect of perturbation, such as 
disorder, anisotropies or external field, is not only neces- 
sary to compare to real materials but also proves to be 
efficient to discriminate between the competing ground 
states - critical QSL are expected to be easily destabi- 
lized while gapped ones should be more robust. 



To confront theories, a significant step was achieved 
with the synthesis [l7| of the first "structurally per- 
fect" QKHA, ZnCu3(OH)6Cl2, called Herbertsmithite. 
The electron spin moments of Cu^"*" ions (S'=l/2) form 
undistorted kagome planes, well separated from each 
other by diamagnetic Zn^"*" triangular planes which en- 
sure the quasi-two-dimensionality of the magnetic net 




FIG. 1. Structure of Herbertsmithite. In this ideal structure, 
Cu and O atoms each occupy one crystallographic site. 



(Fig. [T|). Despite a sizable antiferromagnetic superex- 
change J = 180(10) K, Herbertsmithite develops no on- 
site magnetization [isl . [l9| and remains in an unfrozen, 
fluctuating state under zero field with short-ranged spin 
correlations down to at least 50 mK (~ J/4000), as 
expected for a liquid state [H, Inelastic neutron 

scattering [3] and nuclear magnetic resonance (NMR) 
2ll |22| investigations consistently point to the absence 
of a spin gap larger than ~ J/200 and evidence a peculiar 
power-law temperature dependence of the spin dynamics, 
a possible hint of quantum criticality. At such a low en- 
ergy scale, small perturbations to the ideal nearest neigh- 
bour Heisenberg model may become relevant and deter- 
mine the properties of the ground state. First, intersite 
Cu^+/Zn^+ mixing defects are known to be present at the 
level of 5 — 10 % [is', which could stabilize a gapless 
valence-bond-glass 24]. More fundamentally, the lack 
of inversion center of the magnetic bonds on the kagome 
lattice allows for Dzyaloshinskii-Moriya (DM) interaction 
of spin-orbit origin. In Herbertsmithite, the DM magni- 
tude is estimated to be - 0.04 J - 0.08 J 0, 0], not 
far from a theoretically expected quantum critical point 
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FIG. 2. ^■^O nuclear spin-lattice relaxation rate, l/Ti, as a 
function of T in Herbertsmithite down to 45 mK. Inset: solid 
line represents the initial slope below Te- 



at ~ O.IJ [27|. 

In this Letter, we report the spin dynamics of Hcrt- 
smithite in a magnetic field down to very low tempera- 
ture, T ~ J/4000, probed by ^'^O NMR, in order to un- 
cover the ultimate ground-state properties and test for 
criticality. The location of oxygen on the superexchange 
path between two adjacent Cu ions in the kagome lattice 
confers on the oxygen nucleus a dominant coupling to the 
magnetic properties of the kagome plane (Fig. [ij. This 
site-selectivity is decisive in the context where the low- 
T thermodynamic properties and the response of less- 
coupled local probes are dominated by interplane defects 
[28, 29^. We uncover a field induced transition to a frozen 
phase for applied fields higher than fioHc = 1.55(25) T 
and draw the H — T phase diagram. 

NMR measurements were performed on the powder 
sample of ^^0-enriched Herbertsmithite used in ref. 
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The spin-echo was obtained with the standard pulse se- 
quence, 7r/2 — r — TT. The NMR spectra were obtained by 
recording the integrated echo intensity, point-by-point, 
while sweeping the external magnetic field with a fixed 
rf frequency (30| • The spin-lattice relaxation rate, 1 /Ti , 
was obtained by saturation-recovery method. The abso- 
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FIG. 3. H-T phase diagram of Herbertsmithite from this 
study (black circles) and from the absence of spin freezing 
at least down to 50 mK in /iSR experiment (black square) 
[igl ]. Tc separates the spin-liquid phase from the spin-solid 
phase characterized by strongly suppressed fluctuations. The 
dashed line is a plot of {H — HcY with fioHc = 1.53 T and 
S = 0.65 (see text). Inset: A versus Tc. 



lute values of Ti were obtained by scaling the recovery 
curves to previous results at 6.8 T (2ll. Isoj. 

In Fig. [21 we display our 1/Ti results as a function 
of temperature down to 45 mK in four different mag- 
netic fields. At fi^H = 6.8 T, in contrast to the known 
modest sublinear T-dependence of l/Ti ^ 7^0.73 ^i^Qve 
1 K [2l[, we find a three orders of magnitude decrease 
of 1/Ti below Tc = 520(30) mK. Both the marked onset 
and the magnitude of the drastic slowing down of the un- 
derlying dynamics indicate a transition from the high-T 
spin-liquid phase to a low-T solid-like phase. This dras- 
tic change of the dynamics is detected not only in the 
vicinity of the peak of the NMR spectrum but also at 
various positions on the wings, which points to a homo- 
geneous transition. The 1/Ti data taken at 2.0, 3.1 and 
12.0 T display a similar drop below a temperature T^. 
which is found to increase with the applied field. The 
corresponding phase diagram is drawn in Fig. [3] By ex- 
trapolation to still lower fields, it seems that Tc vanishes 
at some small but finite critical field He- This is in line 
with the absence of freezing in zero field evidenced with 
/iSR [IB] • Fitting the data to a phenomenological expres- 
sion Tc (X {H — HcY yields a quantum critical point at 
HoHc = 1.55(25) T and (5 = 0.6(1). Note that removing 
the 12 T data point, which may be too far away from 
the critical regime, improves significantly the quality of 
the fit (see dashed line in Fig. The Zeeman energy 
Mb^c ~ J/180, where fiB is the Bohr magneton, is quite 
small on the scale of J which underlines the fragility of 
the zero-field spin-liquid state of Herbertsmithite. 

The inset of Fig. [5] shows an Arrhenius plot of 1 /Ti 
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versus 1/T. The linear behavior below T^, evident at 
6.8 T and 12.0 T, strongly supports the existence of a 
gap in the spin dynamics. This gap in the spin-solid 
phase likely reveals some magnetic stiffness arising from 
the DM anisotropic interaction. For T < 200 mK, the 
relaxation originates from a residual contribution which 
is marginal at high fields but sizeable at ij^qH = 2.0 and 
3.1 T. This residual relaxation may be attributed to the 
interplane Cu defects which behave like nearly free S — 
1/2 spins with a small antiferromagnetic coupling of ~ 
1 K [23[ . Their contribution is expected to decrease as the 
field increases, i.e. as their magnetization gets saturated. 
From the initial slope of the l/Ti curves just below Tc 
(Fig. [21 inset), discarding the lowest T data points, we 
extract the gap A. The inset of Fig. [3] displays the values 
obtained for A and their comparison to Tc . Within error 
bars, A increases linearly with Tc with the slope a — 
A/keTc = 2.3(3) (H. 

We now turn to the NMR spectra which reflect the 
distribution of local magnetization around the probe nu- 
clei. For > 3 T and for T > T^, we find that the 
spectrum width is both field and temperature indepen- 
dent (see Fig. 2]). This is in agreement with an earlier 
study at fj,QH ^ 7 T [2lj which demonstrated a size- 
able broadening of the spectra with 1/T at high tem- 
perature and a saturation below 5 K. It suggests that 
the width is related to the interplane defects above Tc. 
The smaller width observed at 2 T and T = 0.2 K is 
in agreement with this picture when the defects are be- 
coming less saturated. At this smallest field, the width 
cannot be measured reliably above T > 0.2 K because 
of inhomogeneous, strong spin-spin relaxation processes. 
Spectra at 12 T, where the longer spin-spin relaxations 
do not affect the lineshape, are displayed in Fig. [5] Below 
Tc, the increase of the width, much like that of an order 
parameter, contrasts sharply with the absence of vari- 
ation observed above Tc. This points clearly to a static 
magnetization which develops in the kagome plane in the 
spin-solid phase. The Gaussian-like lineshape implies a 
random static phase whereas, in the case of an ordered 
phase, the well-defined internal field at the O site would 
yield a rectangular shape for a powder averaged spec- 
trum [sij. Given the hyperfine constant 3.5 T/^b, and 
independently from any model, we can extract a typi- 
cal moment value ficn 0.1/xb, an order of magnitude 
smaller than the full moment for a spin 1/2. This small 
value likely refiects a renormalization of the moment by 
strong quantum fluctuations. The absence of any clear 
broadening below Tc for the lower flelds suggests even 
weaker Cu moments, in line with an enhanced quantum 
reduction in the vicinity of the critical point He. 

We now aim to discuss which model could best account 
for a weak-field induced transition in the context of quan- 
tum criticality. The algebraic (critical) spin liquid model 
[l^ captures both the absence of a gap in zero field and 
the instability under an applied field. A transition to 




FIG. 4. Linewidth (full-width at half-maximum) of the ^^O 
NMR spectra vs. T for the different fields investigated here, 
after correction for inhomogeneous spin-spin relaxation effect. 
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FIG. 5. ^^O NMR spectra normalized to their maximum in- 
tensity in ~12 T at different T. The vertical arrow shows the 
location at which l/Ti measurements were performed. The 
dashed rectangle represents the expected spectral shape for 
120° long-range-ordered (LRO) phase with O.lfiB static mo- 
ments. 



an XY long-range ordered phase and an accompanying 
on-site magnetization M oc was predicted to occur 
at Tc (X H [3]. This is not consistent with our data 
which suggest a finite critical field. Also, the presence of 
a sizeable DM interaction in Herbertsmithite should in- 
duce long-range order already in zero field in this model 

A different approach using exact diagonalization 27] or 
Schwinger boson techniques [H, 34| consists, on the con- 
trary, of including the DM interaction in the Heisenberg 
model as a starting point. Although the QKHA model 
may be gapped [7], the DM perturbation might natu- 
rally account for the gapless zero-field phase of Herbert- 
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smithite by admixing singlet and triplet states. Also in 
this framework, Herbertsmithite lies close to a quantum 
critical point that is in fact sustained by DM anisotropy. 
A quantum critical region is expected to open at finite 
temperatures within a finite range of a control parameter 
around its critical value (s^l , which could in this case be 
the magnetic field. From the point of view of symmetry- 
breaking, the effects of the DM and Zeeman terms may in 
essence not be very different. Such an interplay between 
these two terms was established for S = 1/2 Heisenberg 
antiferromagnetic chains, a canonical example of a QSL 
state in one dimension 36| . Then, Herbertsmithite could 
be driven from a QSL to a spin-solid phase by a moder- 
ate field. He ~ Dc — D. In this framework, the smallness 
of the moment expresses the weakness of the order pa- 
rameter in the vicinity of the quantum critical point [27| . 
Moreover, the quantum critical region is characterized 
by the absence of energy scale other than temperature 
which can explain the power-law behavior of spin corre- 
lations 18, 23-22[. It is then interesting to compare the 
results with other kagome compounds of difi^erent D/J, 
e.g. vesignieite, with possibly D/J > 0.1, shows an un- 
usual heterogeneous ground state where spin liquid and 
frozen moments coexist 37 , [s^ . 

Our results share certain features with 
the field-induced transition recently reported 
for the organic triangular quantum antifer- 
romagnets k-(BEDT - TTF)2Cu2(CN)3 and 
EtMe3Sb[Pd(dmit)2]2 which lie close to Mott tran- 
sition and show a QSL ground state [s^ H^]. As 
in Herbertsmithite, small magnetic fields induce the 
freezing of the zero-field QSL ground state toward a 
spin-solid phase with strongly suppressed moments. 
In NMR measurements, no critical divergence in 1/Ti 
nor critical broadening in the spectra are detected 
43-42 1, similarly to Herbertsmithite as shown here. 
This similarity is quite intriguing since the two lattices, 
triangular and kagome, have been believed for long to 
yield very difi^erent ground-states and excitation spectra 
according to numerical studies 0, [^, [1] . 

In conclusion, the QSL ground state of Herbertsmithite 
appears quite fragile against external magnetic field. 
However, in Herbertsmithite, the observed spin freezing 
likely results from the combination of both a finite ap- 
plied field and sizeable DM interaction which eventually 
suggests that the QSL ground state of pure QKHA it- 
self is somewhat robust against perturbations. It opens 
up the possibility for a spin gap in the pure QKHA in 
accordance with most recent calculations 0. The field 
induced frozen phase that we evidenced here calls for 
further investigation. Experiments in higher magnetic 
fields that saturate the interplanar defect contributions 
and induce larger moments within the kagome plane are 
desirable. Our finding along with the results reported 
for the quantum triangular antiferromagnets, prompts a 
comparative study of the growing class of QSL materials 



to understand their generic behavior in a broader context 
of frustration and quantum criticality. 
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